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Abstract 

In this note we investigate the last passage percolation model in the presence of macroscopic inho- 

mogeneity. We analyze how this affects the scaling limit of the passage time, leading to a variational 

problem that provides an ODE for the deterministic limiting shape of the maximal path. We obtain 

a sufficient analytical condition for uniqueness of the solution for the variational problem. Conse- 

,S^ • quences for the totally asymmetric simple exclusion process are discussed. 

-)— » 
a . 
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1 The model and results 
> 

O. 

The last passage percolation process has been widely studied over the last few years [U [2l [3l [7l [H 
OO ■ [51 [TU]. There are several equivalent physical interpretations for the model. Examples include zero- 

temperature directed polymer in a random environment, a certain growth process, queuing theory, a 
^^ ' randomly increasing Young diagram and random partitions [51 1121 \TJ\ . By a simple coupling argument, 

C^ , results obtained for last passage percolation have their duals for the totally asymmetric simple exclusion 

process and thus the former model is often useful for the study of the later one [H O [Til HSl [HI [S] . It 
is for the two-dimensional case that the most explicit results and estimates are known. In particular for 
geometric or exponential distributions, to which an exact solution was given by [H]. 

In this note we study last passage percolation with exponentially distributed passage times in the presence 
of macroscopic inhomogeneity. We begin by restating known results with a different point of view: instead 
of taking the limit of a large rectangle on the usual lattice we consider the limit of a fine lattice on some 
fixed rectangle, which is equivalent but resembles hydrodynamics. A continuous function a is defined on 
the macroscopic rectangle and locally modifies the parameter of the process. The problem of studying 
the random microscopic path of maximal passage time leads to the variational problem of finding a 
deterministic macroscopic curve maximizing a certain functional. We shall see that the rescaled passage 
time indeed converges to that given by the variational problem, and give sufficient analytical conditions 
for convergence of the maximal path's shape to a deterministic curve. 

The viewpoint adopted here has immediate implications for the totally asymmetric simple exclusion 
process. On the scaling limit we can describe the behavior of the total current through the origin up to 
a given time when the jump rate has macroscopic fluctuations in space (as well as in blocks of particles, 
or even both). In particular for spatial inhomogeneity it is easy to see that the instantaneous current is 
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non increasing in time. The analysis of the time taken for a given amount of particles to cross the origin 
and of the corresponding path that gives that passage time helps understanding the bottleneck, that is, 
which chain of events was responsible for that delay. The relation with TASEP will be discussed further 
in Sectional 

The last passage percolation problem may be formulated as follows. Given the origin and a point {l,b) € 
R^, l&l < I, consider the rectangle Q = {{x,y) : ^ x ^ l,\y\ ^ x,\b — y\ ^ I — x], that is, the rectangle 
like displayed in Figure [T] having (0,0) and {l,b) as vertices. For N £N, take the grid Sn = Q n -^Z^, 
where Z^ = {(x, j/) G Z^ : x + y E 2Z}. Let ^n denote a field {£.p)p£Sn: whose coordinates are i.i.d. 
distributed as -^exp(l), a > 0. One should think of £,p as a reward standing at the site p. Let Hn 
denote the set of oriented paths {po, ■ ■ ■ ,Pk) crossing Sn, i.e., po = 0, Pi & Sn, Pi — Pi-i — ;^(l,il) 
for i — 1, . . . ,k and pk is the rightmost point of Sn ■ If tt G Hn , we define ^at • tt as the sum of the 
random variables ^p for p & it, i.e., the total value of the rewards along that path. Let tt^ denote the (a.s. 
unique) path that maximizes £,n ■ i" and let G{£,n) = maxTreiiN (^a' ' t^) = (,n ■ ^^ ■ It is known [8J that 
the maximal value G{S,n) approaches aZ7(j), where 7(w) = 1 + \/l — w'^, as N increases. Moreover, 
the probability of deviating e from this limit decays exponentially fast in N. A consequence of this fact 
is that the maximal path tt^ approaches the straight line connecting the origin and {I, b) in the || • ||sup- 

We now describe how the macroscopic inhomogeneity is introduced. Let a be a nonnegative, continuous 
function defined on Q and for each A^ G N we take {^p)peSN as independent random variables, each one 
distributed as ^^^ exp(l). We are interested in understanding the value of G{S,n) and the shape of tt^ 
for large values of N. 

Let us give some simple heuristic arguments. Take some y(-) G X, where 

A = {y : [0, Z] — > R I y(0) = 0, y{l) — b,y is Lipschitz with constant l}. 

We first look at the macroscopic limit of max^(^jv • ti") restricted to paths tt that stay "pinned" to the 
curve y. Take Ax ^ I, consider the set 11^ C XIat of paths that pass by [x, y{x)) for all x — nAx, and 
let G'{^n) = maxTren' {^n ■ tt). Since a is nearly constant on a small neighborhood of {x,y), we expect 
the contribution to G' obtained between x and x + Ax to be given by a{x,y)^(-^^Ax, in accordance 
with the constant case discussed above. This motivates the definition of the functional 

Q{y)^ / OL{x,y{x))-i{y)dx, y € X, 
Jo 

where y = dy/dx. It is reasonable to expect that G'{^n) will approach G{y) when A^Aa; — > oo sufficiently 

fast and Ax — > 0. Informally we say that that the maximal path n, restricted to be pinned to the curve 

y, will asymptotically catch a total of G{y) in rewards. Now drop this restriction, i.e., take the maximal 

path among all possible paths instead of being pinned to some specific curve y. We guess G{^n) will be 

given by maximizing over y, that is, 

G{^n) ~ sup g{y). 
vex 

Our first result confirms the above argument. 

Theorem 1. Let a be a continuous function on the rectangular domain Q. Consider the inhomogeneous 
last passage percolation problem as described above and take 

Q* = snpg{y). 

yex 

Then G{£,n) -^ G* a.s. Moreover, for any S > there are c, C > such that 

P{\G{U)-G*\>5)<Ce-''' 
holds for all N en. 



It is clear that the variational problem defining Q* is crucial for the understanding of the process. Unique- 
ness of its maximizer is required in order to establish convergence of the maximal paths. Although one 
can easily build a function a for which two distinct curves maximize Q*, the authors believe that the set 
of a's for which y* is unique is generic in the C'^{Q) topology. The next result concerns the shape of the 
maximal path. 

Theorem 2. For a continuous there is y* & X such that Q* = Q{y*). If such y* is unique, the random 
maximal path tt^ approaches the curve y* for N large. That is, \\tt^ — y*||sup — > a.s., and, for any 
(5 > there are c, C > such that 

holds for allN eN. 

Suppose also that a has continuous derivatives ax and ay. Then the Euler- Lagrange equation associated 
with the variational problem of Q* 

2/ = -^ K(i - y')^^' + Ky + «y)(i - f)] n ^ 

y{0) = 0, y{l)^b ^'> 

has at least one solution yo. If in addition a is such that Q is strictly concave on X, the solution ya of |IP 
is unique and it is also the unique maximizer of Q . 

When a is constant it is easy to see that the maximizer must be a straight line, some other nontrivial 
examples of uniqueness are given by the condition we discuss below. Of course uniqueness holds whenever 
Q is strictly concave. When a is smooth, a sufficient condition for strict concavity of Q is 



Yyy < and - aayy > -{ayf (2) 



for all (cc, y) £ Q. 

In Section [2] we prove Theorem [T] with the aid of large deviation estimates already known to hold when 
a is constant. Section [3] contains the proofs of Theorem [2] and the sufficiency of the condition ([2]). In 
Section [4] we discuss some implications of these results for the TASEP. 

The authors believe that it is possible to prove similar results for the case of geometric distribution of 
passage times and for the case of non-homogeneous poissonian points with the same arguments presented 
in this note. 
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2 Convergence of the passage time 

In this section we sketch the proof of Theorem [1] The proof consists on approximating a by a constant 
on smaller regions of the domain Q and then applying the large deviation principle known to hold on 
these regions. 

Since a will be approximated by other functions, we write Qa{y) and ^q • tt to specify which a is being 
considered. Besides, processes with different a's can be constructed on the same probabihty space in a 
way that S,a ^ ^a whenever a^ a. 

We first prove that P(G{£,n) < G* - e) < Ce'"^ . Given e > 0, divide Q as in Figure [TJ the resulting 
squares being fine enough so that the value of a does not change more than ^ inside each of them. Let 



a be constant on each square and given by the infimuni of a on that square. In this case a— -^ ^ a ^ a. 
Take y E X satisfying Giv) > G* — e/3. Since 7 ^ 2 we have 



Gaiy) ^ Gciy) -^^G* 



2e 



Let Qi, . . . ,Qk be the rectangles that contain the curve y between its intersections with the grid as in 
Figure [TJ GaiQj) stands for the supremum of G restricted to Qj (the functions and the integration are 
both restricted to this set, Qj is a rectangle and the functions connect both extremes). 





Figure 1: On the left, the definition of the sets Qj of a given y. On the right, the k strips of Q and a 
given random path tt . 



Take 11^ C n^v as the set of paths that pass through the left and the rightmost points of every Qj . Given 
TT S n^ write ttj for the restriction n to Qj. If (^ • n^) < G* — e, the same must hold for every n in 11'^ 



N 



(recall that n is the maximizer of S, ■ tt). This implies that for some j, ^5 • ttj ^ ^^ • ttj < G^iQj) — e/3fc 
for every tt in 11^. But by [14j the probability of such event decays exponentially. So, 



P[Ca -n^ <G* -e] s; fcmaxQe""^^ < Ce 



-cN 



for suitable C and c. We point out that in fact P\_^a ■ tt^ < G* ~ e] ^ Ce '^^ , since the large deviation 
principle proved by [Til is on N'^ and not N. 

Finally let us show that P[G{£_n) > G*+e) < Ce""'^ . Given e > 0, take S > such that \a{pi) - a{p2)\ < 
e/8l whenever \pi — P2I < 2(5. Divide Q in fc strips of width r = l/k as in Figure[Tl with k such that r < 5, 
and for each path vr G IIjv write tt^ for the restriction of tt to the j-tli strip (note that S, ■ tt — ^,S, ■ ttj). 

For a given tt e 11, take pi = (xi, j/i), . . . ,pfe^i — (xfc_i, yk-i) as the points where the trajectory of tt 
intersects the boundary of each strip and let yk = b. Note that pj needs not to be on the lattice. 

We will use the following consequence of Theorem 1.6 in 8 . If a ^ M is constant, there are C,c > 0, 
depending only on M and e, and such that, for any choice of j, yj and yj+i, 



P [37rj, ^ • TTj > ar-i{{yj+i - yj)/r) + e] < Ce 



-cN 



Denote by 5*[pi,P2] the supremum of Ga taken over the 1-Lipschitz functions that connect pi with p2, 
the integration done on the appropriate domain. Now 



fc-i 



fc-i 

E 



Gl >^Gl[Pj,Pj+i\ > Yl H^i) + ^ Ti{{y]+i 



%)A) - 2' 



The first inequality is obvious. The second one follows by the choice of S, since it implies that for any 
1-Lipschitz function y connecting pj and Pj+i we have a(x, y{x)) ^ C({pj) — -^ and 7^2. 

Therefore, if [^a ■ tt^] > t/* + e, there must be j such that £_a ■ ttj > £,a • t^j ^ drj(^{yj+i — yj)/r) + ^, 
where a = a{pj) + ^ satisfies a > a{p) for all p G ttj. Considering all possible choices of j, yj, yj+i we 
get 

P[Ga{0 > a* + e] s; IMkN'^Ce-''^ < Ce"^^, 

concluding the proof of Theorem [TJ 

3 Limiting shape of the maximal path 

In this section we present the proofs of Theorem [2] and the sufficiency of ([2]) for uniqueness. 

We first address the question of existence of a maximizer y* . By Ascoli-Arzela, X is a compact space 
with the uniform topology. However, Q is not a continuous functional in that topology (notice that Q{y) 
is defined in terms of the derivative of y, which exists a.e. in [0, 1] for y S X). Nevertheless the proof will 
follow by compactness, since the Q is upper semi-continuous, as we prove now. 

We shall show that Q = inf„i Q,n, where the "Riemann sums" 

are clearly continuous in the uniform topology if a is continuous on Q. Since 7 is a concave function, 
it follows from Jensen's inequality that Q ^ Qm for every tti, so it remains to show that Qmiy) ~^ Giy) 
for all y. Let y £ X and take a sequence of y„ £ C^ such that ||j/ — yn||oo — *■ and \\y — ynllii — > as 
71 — !■ 00. Write 

1^(2;) - 5m(y)| =^ \Q{y) ~ S(2/«)l + \Q{yn) - Qm{yn)\ + \Qm.{yn) - Qm.{y)\- 

The first term vanishes as 7i — > cx3 by the bounded convergence theorem. The second term vanishes as 
771 ^ 00 for n fixed, since for smooth y the numbers Gm{y) indeed correspond to Riemann sums of a 
continuous integrand along the path x 1— > (x, y, y). It thus suffices to show that the last term goes to zero 
uniformly in ?7i as n ^ 00. Since Ijiw) — ^{w')\ ^ C-\/|w — w'\ for all w, w' , 



\Qmiyn)~Gmiy)\ ^Ca^ — i / \y{lz/m) -ynilz/m)\dz] ^ Cay/l\\y-y: 

771 \ Ji_l / 



where the second inequality is Jensen's inequality for the square root and a denotes supg |a|. This proves 
that G is upper semi-continuous, therefore there is a y* attaining supj^g^ Giy)- 

We now discuss the convergence of tt^ to y* when this is the unique maximizer. It follows from uniqueness 
of y* and upper semi-continuity that 

for any S > 0, 3e > such that y £ X, \\y — 7/*||sup > 5 imphes Q{y) < Q* — e. 

This, together with the exponentially fast convergence in probability given by Theorem [l] implies the 
desired result. 

Suppose a has continuous derivatives and consider the Euler-Lagrange problem ([T]). For a given C^ 
function y £ X with \y\ < 1, it is easy to see that ([T]) holds if and only if 

§:g{y + th)^ =0 

^'' lt=o 



for every smooth perturbation h that vanishes on the extremes of the interval. Rewrite ^ as 

^ = ~i K(l " w'^f/'^ + {a^w + ay){l - w^)] 
y = w 

and consider the initial value problem given by y(0) = 0, w(0) = Wq and ([3]). Existence, uniqueness and 
continuous dependence of solutions on the initial condition wq follow from basic ODE theory. Now notice 
that for Wo = ±1 we have w{x) — ±1 Vx and y{l) = ±1. So, for given b £ (—1,1), by the intermediate 
value theorem there is wq £ (—1,1) for which the corresponding solution j/o satisfies yo{l) — b. Since 
{w = ±1} is an invariant set, ya must satisfy |yo(a;)| < 1 V a:. Therefore there is at least one solution 
2/0 e ^ of H]). 

Finally we show that when Q is strictly concave, the solution yo is unique and it is also the unique 
maximizer oiQ. We simply adapt a standard method for showing that a critical point of concave functional 
must be the unique maximizer. Let yQ £ X he a solution of ([1]). Suppose there were j/i £ X\{yo} with 
Giyi) ^ Giyo)- By the strict concavity, there is y2 S -'^ such that ^(2/2) > Giyo)- Approximate 2/2 by 
some smooth y^ £ X such that still Giys) > Giyo)- The concavity of G imphes ^^(2/0 + tfT')\t=a > 0' 
where h — y^ — yo, a contradiction. Therefore 2/0 is the unique maximizer of t/ in X. As 2/0 was taken as 
any solution of the contour value problem, this must be the unique solution. 

This finishes the proof of Theorem [2l We end this section showing that ([2|) is sufBcient for G to be strictly 
concave. 

Let xo £ (0,1) be fixed and define a{y) = a{xo,y), z{y,w) = a{y)j{w). What we shall actually prove is 
that the function z is strictly concave, which is much stronger than the functional G having this property. 
We study the eigenvalues of the Hessian H^ of z to have conditions under which z is a concave function. 

det(iJ^ - A) = A^ - {aj" + a"j)X + [aa"jj" - {a'j'f]. 



Both eigenvalues are negative if and only if s < and p > 0. Since 

(7')' 1 

it is enough to require ^ to have the strict concavity of z and consequently of Q. 

4 Applications to TASEP 

The last passage percolation model can be coupled with the totally asymmetric simple exclusion process 
with the initial condition that the sites are occupied iff they lie to the left of the origin. See [15] for a 
description of such correspondence. The current through the origin for the TASEP can be understood 
by considering last passage percolation between (0, 0) and {I, 0), since this gives the rescaled time needed 
for lN/2 particles to cross the origin. At the scaling limit this converges to G*[l]- 

The model considered in the previous sections can describe, for example, the TASEP with the jump rates 
depending smoothly on the position of each site, when we take a depending on y. In the hydrodynamic 
limit, the instantaneous current through the origin will be a decreasing function of the time. Intuitively, 
as the time passes, the process finds more bottlenecks that were not important before, and at some point 
such barriers start being determinant for the passage time, at least until the system meets even narrower 
bottlenecks. 

The above fact follows easily from the variational formulation of G- Consider for instance the graph of 
a{y) having two peaks as in Figure [21 For short times the system dynamics does not feel the existence 





Figure 2: on the left an example of a{y) having two strong peaks; on the right the maximizers ?/* for I in 
three different regions. 



of such regions of large average jump time, so Q*[l] (i.e., the time needed for Nl/2 particles to cross the 
origin, at the scaling limit) grows linearly with I as in the homogeneous case. For larger times the first 
peak will become important: there will be a traffic jam before this point and the rescaled time G*[l] will 
increase roughly linearly with /, but at bigger rate, given by the value of a at this peak. Now for even 
larger times again the system finds a harder difficulty to overcome, which is to wait for particles to cross 
the stronger bottleneck given by the second peak, leading to a low density of particles flowing through 
the origin, so from this time on Q* [I] increases at an even higher rate, again given by the value of a at 
the new peak. The curves shown in Figure [2] give the maximizers y* of 5*[^], for I in these three regions, 
illustrating the decreasing of instantaneous current. 

For a proof of this phenomenon notice that the instantaneous current being non increasing is equivalent 
to saying that the amount of time needed for a certain flux to pass by the origin is a convex function. 
Now the following inequality implies that / t—f G*[l] is convex. Given Zq > 0, 

g*[l]^g*[lo] + 2a{l-lo), 

where a = snp,Q iTa(y{x)) — Q;(y(xo)) and G{y) = ^*[^o]- To see why this inequality hods for I > Iq 
we translate the curve y hy I — Iq after xq and fill the gap with a constant line; for the resulting curve y 
we have g*[l] > g{y) = Q{y) + 2a{l — Iq) — G*[lo] + 2a{l — Iq)- For I < Iq we find a point xi such that 
yi^i) = yi^i +lo~l) and concatenate y\[o,xi] with y\ixi+io-i,io]- 
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